Variation of the glass transition temperature with rigidity and chemical composition 
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The effects of rigidity and chemical composition in the variation of the glass transition temper- 
ature are obtained by using the Lindemann criteria, that relates melting temperature with atomic 
vibrations. Using this criteria and that floppy modes at low frequencies enhance in a consider- 
able way the average cuadratic displacement, we show that the consequence is a modified glass 
transition temperature. This approach allows to obtain in a simple way the empirically modified 
Gibbs-DiMarzio law, which has been widely used in chalcogenide glasses to fit the changes in the 
glass transition temperature with the chemical composition . The method predicts that the constant 
that appears in the law depends upon the ratio of two characteristic frequencies (or temperatures). 
Then, the constant for the Se\- x - y Ge x As y glass is estimated by using the experimental density of 
vibrational states, and the result shows a very good agreement with the experimental fit from glass 
transition temperature variation. 

PACS numbers: 64.70.Pf, 64.60.-i, 05.70.-a, 05.65.+b 



Glasses are solids that do not have long range order, 
and which are usually produced by fast cooling of a liquid 
melt. In spite of the well known importance of these ma- 
terials in human history, the physics of glass formation 
is still an open and puzzling problem p|. The reason 
behind this situation is that glass formation is mainly 
a non-equilibrium process, and thus even some impor- 
tant questions for glass technology are not well under- 
stood, as the origin of the non-exponential relaxation 
laws [2J or the ability of certain materials to reach the 
glassy state Q. The first approach to the problem was 
made by Kauzmann Q, who pointed out that there 
is an underlying thermodynamical phase transition be- 
hind glass formation due to an entropy crisis. Although 
the Kauzmann approach was very successful, it was un- 
able to answer many questions, specially in a quantita- 
tive way. Later on, other approaches were proposed p|: 
phenomenological models like the Gibbs-Dimarzio, the- 
ories based in supercooled liquids as for example mode 
coupling and energy landscape formalism, or the use of 
extensive computer simulations 0. However, most of 
these theories have difficulties to explain one of the most 
simple and interesting questions: how the glass transi- 
tion temperature (T g ) depends on chemical composition. 
As discovered by the Phoenicians more than 2,000 years 
ago, T g can be dramatically lowered or raised by adding 
few impurities. Another interesting property of glasses 
is the behavior of viscosity, which is usually referred as 
fragility [|| . Fragility is related to the glass forming ten- 
dency, since a non-fragile glass former do not require a 
high speed of cooling to be produced. The fragility of the 
glass can be changed from strong to fragile Chalco- 
genide glasses (formed with elements from the VI col- 
umn doped with impurities) are very important to un- 
derstand these effects. In general, the change of T g with 
the chemical composition has been observed to follow an 
empirically modified Gibbs-DiMarzio law 6] 7] 8] 9] . It is 
worthwhile mentioning that the original Gibbs-DiMarzio 



law was developed for polymers |4[, and contains some 
special assumptions that have been very elusive to mea- 
sure, like the size of cooperative rearranging regions Q. 
Furthermore, the use of this law for chalcogenide glasses 
is based in the hypotesis that these glasses nearly behave 
as polymers. A completely different approach is to use 
stochastic matrices [l(J, |l2j . The problem with this 
theory is that one needs to assume certain conditions to 
define the glass transition which are not in the form of 
usual definitions [ToT ]. 

The rigidity theory (RT), introduced by Phillips [l3| 
and refined by Thorpe |l4j . was another important step 
to understand glass formation. By considering the cova- 
lent bonding as a mechanical constraint, the ease of glass 
formation is related with the ratio between available de- 
grees of freedom and the number of constraints. If the 
number of constraints is lower than the degrees of free- 
dom, there are zero frequency vibrational modes called 
floppy 15]. The resulting network is under-constrained. 
A transition occurs when a disordered lattice becomes 
rigid. Glasses that are rigid at a certain chemical compo- 
sition are easier to form, and many features of this transi- 
tion have been experimentally observed [g] 16]. Even for 
simple systems like hard-disks polymer melts [l8j . 
and colloids [19j, rigidity plays an important role . For 
more complex systems like proteins, RT is a very powerful 
tool to understand folding and long-time scale motions 
po|. Recently, it has been discovered that at the rigid- 
ity transition, there is a window of thermodynamical re- 
versibility 8] 21] 22] 23], explained as a self-organization 
of the stress in the network [24| 25]. Although RT ex- 
plains a lot of the phenomenology of chalcogenide glasses, 
the introduction of thermodynamics into RT has been 
elusive 18]. In previous papers, we approached some of 
these problems by using a phenomenological free energy 
|26|. and then we made computer simulations p7||28[. 
Yet, many questions remain open, as for example the 
variation of T g with the chemical composition of the glass. 
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In this letter, we will show that by using the RT and the 
old Lindcmann criteria that relates melting with the size 
of atomic vibrations (which has been proved to be valid 
for glasses H^ll^), one can explain the variation of T g 
with chemical composition, and obtain in a simple and 
elegant way the modified empirical Gibbs-DiMarzio law. 
Furthermore, it provides a clear path to understand glass 
transition trends with the connectivity of the network. 

In RT, the ability for making a glass is optimized when 
the number of freedom degrees, 3N, where N is the num- 
ber of particles, is equal to the number of mechanical 
constraints (N c ) that are given by bond bending and 
stretching forces. The parameter / = (37V — N C )/3N 
gives the fraction of cyclic variables of the Hamiltonian. 
A variation of these coordinates do not change the energy 
and thus / also corresponds to the fraction of vibrational 
modes with zero frequency, called floppy modes |14| , with 
respect to the total number of vibrational modes. The 
parameter / is in fact a function of the average chemi- 
cal coordination of the glass, and can be calculated using 
the pebble game algorithm [15| , but usually it can be es- 
timated with a mean-field procedure, known as Maxwell 
counting |3lj . 

3 3N 6 w 

where (r) is the average coordination number, defined as, 



(r) = rx r 

r=0 

and r is the coordination number of each atomic specie. 
x r is the occurrence of each type of atom [3lJ . A rigid- 
ity transition occurs when / = 0, and the system passes 
from a floppy network to a rigid one. Glasses are rigid if 
there are more constraints than degrees of freedom, and 
then the lattice is overconstrained. In 3D, the rigidity 
transition leads to the critical value (r c ) = 2.4 if all an- 
gular constraints are considered |3lj . In principle, since 
floppy modes have zero frequency; they do not contribute 
to the elastic energy, and the specific heat of the glass de- 
sends on /, a result that violates the Dulong-Petit law 
26(, and which is not observed in the experimental data 
g . In fact, floppy modes do not have a perfect zero 
frequency, i.e., in real glasses they are shifted by resid- 
ual forces, like the Van der Waals interaction [Tj]. This 
is confirmed by neutron scattering experiments, where it 
has been shown that floppy modes in the prot otyp ical 
compound Sei- x - y Ge x As y are blue-shifted 0|22] 33], 
forming a peak around bmeV. The corresponding tem- 
perature (6/) where these modes are frozen, can be 
estimated from the energy required to excite modes of 
5meV, that gives 0/ ~ 54° K. For low temperatures, is 
clear that floppy modes are important, as confirmed by 
the giant-softening of the 1195'n Lamb-Mossbauer factor 
in (Geo.ggSYio.oOzS'ei-z glasses [l^. But in principle, 
one can argue that floppy modes are not so important 
for the physics of the glass at high temperatures, since 



all 37V/ floppy oscillators are already excited, leading to 
the Dulong-Petit law. As we will show next, this is not 
the case, and in fact floppy modes are essential in deter- 
mining T g . To show this, we use the Lindemann criteria 
(1910), which has been originally devised to understand 
crystal melting 34]. Such criteria establishes that melt- 
ing occurs when the mean atomic square displacement 
(u 2 (T)) is around 10% of the atomic spacing a. There 
are many experimental evidences that this criteria is ap- 
plicable to glasses |3CL as for example happens with Se 
and Se - Ge systems @- The value of (u 2 (T)) can be 
calculated from the density of vibrational states g{u>). At 
high temperatures, 



(u 2 (T)) 



mJV 



(1) 



where k is the Boltzmann constant, T the tempera- 
ture, lu the frequency and m the mass. At T g , the Lin- 
demann criteria applied to glasses [2^] establishes that 
(u 2 {T g )) m (u 2 {T m )) « 0.1a 2 where T m is the melting 
temperature. This shows the fundamental importance 
of floppy modes to determine T g due to the enhance- 
ment of low frequencies produced by the term uj~ 2 in 
eq.JTJ, which leads to an increasing (it 2 ) .Let us made 
a model for the most prototypical chalcogenide glass: 
Sei^ x - y Ge x As y , since this compound allows to obtain 
the same (r) with different chemical compositions. The 
three chemical elements of this glass have nearly the same 
mass; thus we can suppose that in eq.JIJ all the atomic 
masses are equal, although in more general cases m must 
be replaced by an averaged mass. The most simple form 
for g{u)) is to use an Einstein like model, with a delta 
function centered around a characteristic floppy peak at 
LOf, with a weight given by the number of floppy modes, 
plus a density of states that has the rest of the spectral 
weight. Such a density of states can be written as, 

g(w) = 3N(1 - f) gi (u) + 3NfS(co - co f ) (2) 

where gi(u>) is a density normalized to one for / = 0. It 
is important to remark that in this model, gi(oj) will be 
the same at all glass compositions below (r) = 2.4, and 
the relative weight on non-floppy modes is taken into ac- 
count by the factor 1 — /. Note that by using a fixed 
gi{io) for all chemical compositions, we are overestimat- 
ing this contribution for flexible glasses, which will give 
only a small correction to the final result, since g{u>)/w 2 
goes to zero at high frequencies [321 . This assumption 
is also supported by the neutron scattering data [33| . 
which shows that g(u>) follows an isocoordination rule for 
Se±- x - y Ge x As y . Using g(ui), (u 2 (T)) turns out to be, 



(u 2 (T)) = 



3kT 




(3) 



where ( -K ) „ is defined as the second inverse moment at 
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the rigidity transition, 

oo 

1 \ f gi(cj) 

R 



Eq.(|3Jl predicts a linear dependence of (m 2 (T)) upon /. 
This result and the idea of an Einstein-like mode for 
floppy modes are supported by the giant-softening of the 
119SYt Lamb-Mossbauer factor in (Geo.ggSno,oi) x Sei- x 
glasses as x linearly decreases to 0, i.e., as the glass 
gets more floppy. The mean square displacement at ab- 
solute zero or the first moment of the vibrational den- 
sity of states contains the information on these floppy 
modes, and its variation x nicely scales with the scat- 
tering strength of the bmeV mode observed in inelastic 
scattering [lq. 

Now we turn our attention to combine these results 
with the Lindemann criteria. When / = 0, we use eq.© 
and the Lindemann criteria to obtain T g at the rigidity 
threshold / = 0, 



Tg(J = 0) 



O.lma 2 



(4) 



For the variation of T g as a function of rigidity, we use the 
previous expression to rewrite eq.© in terms of T g (f = 
0), and then we apply again the Lindemann criteria to 
obtain, 



Tg(f) = 



Tg(f = 0) 

1 + af 



(5) 



where T g (f) is the glass transition temperature when a 
fraction / of floppy modes is present. The parameter a 
is defined as, 



( J 2 )j 



- 1 = 



Or 



1- 



and depends upon the ratio of two characteristic frequen- 
cies, since we define the following frequency and temper- 

I/O 

ature at the rigidity transition, Qr/H = ujr = (jjjr} fl 

Equation (JSJ predicts a decreasing glass temperature 
with the number of floppy modes. Furthermore, by using 
the Maxwell mean field to express / in terms of < r >, 
eq.© can be transformed into, 



T g ((r)) = 



T g ((r) = 2.0) 



1 - /3 «r> - 2.0) ' 
where the constant (3 is given by, 

5 a 



(6) 



2a + 6' 



The value of T g ((r)) obtained in eq.© is exactly in the 
same form of the empirically modified Gibbs-DiMarzio 



equation that has been used as a fitting law for many 
chalcogenide glasses, where (3 is a constant taken from 
the experimental data. Our approach shows that this 
constant depends upon the ratio of two characteristic fre- 
quencies or temperatures. It is worthwhile mentioning 
that the constant experimental value of (3 in the interval 
2.0 < (r) < 2.4 gives an extra support to the assumption 
made in eq.JSJl. For (r) > 2.4, the isocoordination rule 
for T g is broken and (3 is no longer constant for different 
chemical compositions, as is also observed with the glass 
transition temperature. Another advantage of the ap- 
proach presented in this letter, is that (3 can be estimated 
from neutron scattering data or Lamb-Mossbauer factor. 
Here we will use the neutron scattering data taken from 
reference [s^. We start by noting that at the rigidity 
transition, f — and from eq.©, gi(u>) = g(w).Thus, 
(u7 ! )r can ^ e obtained from the experimental data by 
calculating the second inverse moment of the normal- 
ized density of states at (r) = 2.4. Following this pro- 
cedure, we get (-^r) R ~ 0.QlQ%&meV~ 2 , with a char- 
acteristic frequency ujr w 7 .0959meV , and temperature 
Q R sa 76°ii'.Using that ojf w 5meV, a has the approxi- 
mate value 1.014, and finally (3 « 0.67. The value from 
fitting the experimental data of the glass transition has 
been obtained by many groups @ @ , and produce the 
value (3 « 0.72 which is in very good agreement with our 
estimation from neutron data. The calculation of (3 can 
be improved by taking into account that the floppy peak 
is not a delta function, but instead it has a finite width, 
which will provide a small correction to the result. How- 
ever, even in this simple first approximation the results 
are encouraging. 



In conclusion, we have explored the effects of floppy 
modes into the glass transition temperature of glasses. 
Although these modes are not at zero frequency and at 
high temperature they do not play a role in the specific 
heat, they are fundamental to determine the glass transi- 
tion temperature since they enhance the square mean dis- 
placement. Then, by using the Lindemann criteria and a 
very simple model of the vibrational density of states, we 
recovered the empirically modified Gibss-DiMarzio equa- 
tion that has been widely used to account for the change 
in the transition temperature of chalcogenide glasses. 
The constant that appear in this law depends upon two 
parameters: the frequency of floppy modes and the sec- 
ond inverse moment of the density of non-floppy modes. 
The results of this letter are in agreement with Phillip's 
original idea that glass forming tendency is enhanced at 
the rigidity transition [35|[3ii, as also thermodynamical 
reversibility 0, since it seems that there is less entropy 
due to a smaller quadratic displacement. 
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